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1 .  Introduction 

The  research  supported  under  this  grant  was  concerned  with  analytical  and 
numerical  investigations  on  the  basic  semiconductor  equations. 
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The  research,  as  outlined  in  the  original  proposal,  was  grouped  into  two  areas, 
an  analysis  of  the  mathematical  structure  of  solutions  via  a  singular  perturba¬ 
tion  approach,  and  the  development  of  numerical  methods  for  the  transient 
problem.  The  analytical  part  was  divided  into  several  subcategories.  First, 
the  transient  problem  was  analyzed  in  the  case  when  the  mobilities  were  only 
functions  of  position.  Then  velocity  saturation  effects  were  included  so  the 
more  complicated  case,  when  the  mobilities  depend  on  the  electric  field,  was 
treated.  This  was  done  for  one -d imensiona 1  problems  first  and  later  on  extended 
to  the  two-dimensional  cav.  n  the  numerical  part  difference  schemes  were 
developed  which  incorpora?  ••  ?  •  ••  information,  obtained  in  the  first  part,  about 
the  structure  of  the  in\<  :  r  itors.  Later  on,  the  stability  properties  of 

these  schemes  were  anal../’  »•  ’irst  discuss  the  general  progress  in  the  area 
before  explaining  thes>-  •  iare  detail. 


2. 


General  Progress 


Several  results  have  been  published  in  the  last  two  years  regarding  the 


existence  of  solutions  and  the  validity  of  the  singular  perturbation  approach  in 
the  steady  state  case.  In  [7  ]  the  existence  of  solutions  to  the  time  dependent 
problem  has  been  shown.  In  [1  J  the  convergence  of  the  solution  to  that  of  the 
steady  state  problem  has  been  proven  in  the  one-dimensional  case  and  convergence 
rates  have  been  derived.  In  6j  this  has  been  extended  to  the  multidimensional 
case.  Since  this  result  is  derived  using  Schauder  arguments,  however,  it  does 


not  provide  convergence  rates.  A  singular  perturbation  approach  of  a  slightly 
different  flavor  has  been  introduced  in  [3],  [4l.  It  is  suited  for  reverse 
biased  devices  and  allows  the  calculation  of  the  boundary  of  the  charge  free 
region  via  the  solution  of  an  obstacle  problem.  In  [2  ]  a  difference  method  for 
the  transient  problem  has  been  developed  which  has  features  simular  to  the  one 
developed  under  this  contract  (i.e. ,  second  order  and  L-stable). 


Achievements  and  Publications 


Analysis: 

In  ^A1  and  [B  ]  the  proposed  asymptotic  analysis  has  been  carried  out  in  the 
case  of  one  space  dimension  and  field  independent  mobilities.  While  in  [A]  a 
symmetry  assumption  had  been  made,  which  simplified  the  analysis,  [B  ]  deals  with 
the  more  practical  case  of  a  strongly  asymmetric  doping  concentration.  In  IC1 
the  mobilities  have  been  made  dependent  on  the  field  in  the  usual  way  such  that 
the  velocities  u(1^»  Uplfy  'saturate'  (approach  a  limit)  for  ( V^j  +  °°  .  This 
causes  some  technical  difficulties  in  the  derivation  of  the  asymptotic  expansions 
and  introduces  an  additional  time  scale,  which  lies  in  between  the  diffusion  time 
and  the  dielectric  relaxation  time.  Numerical  experiments  have  confirmed  the 
presence  of  this  timescale.  In  [D ]  this  analysis  has  been  extended  to  the  two- 
dimensional  model  of  an  MOS  transistor.  Existence  proofs  for  the  solutions  of 
the  reduced  problem  have  been  given  in  the  case  of  reverse  biased  or  very 
moderately  forward  biased  devices.  It  should  be  pointed  out  that  in  the  case  of 
field  dependent  mobilities  the  layer  equations  fail  to  have  a  solution  if  the 
currents  become  too  large. 

Numerical  Methods: 

As  a  first  step,  in  order  to  justify  the  use  of  standard  methods  like  back¬ 
ward  differences  (BDF ) ,  the  linearized  eigenvalue  problem  has  been  analyzed  in 
(E).  Extending  previous  results  by  Mock,  it  has  been  shown  that  eigenvalues  can 
have  positive  real  parts  away  from  steady  states.  These  real  parts  are  bounded 
from  above,  however,  by  moderate  quantities  proportional  to  the  size  of  the 
currents.  Using  this,  one  immediately  can  derive  stepsize  criteria  for  BDF 
methods  which  guarantee  linearized  Von  Neumann  stability.  In  IF!  and  [C j  a  time 


differencing  scheme  has  been  introduced  and  analyzed  which  exploits  the  special 
structure  of  the  stiffness  of  the  basic  semiconductor  equations.  Similar  to  the 
method  developed  in  2  j  it  is  L-stah).e  and  of  second  order  as  long  as  the  time- 
steps  are  much  larger  than  the  dielectric  relaxation  time.  However,  in  contrast 
with  the  method  in  2],  it  requires  the  solution  of  only  one  nonlinear  system  per 
timestep.  The  task  of  performing  a  nonlinear  stability  and  convergence  analysis 
of  this  method  has  proven  too  difficult  so  far.  However,  using  the  analysis  of 
E  a  linearized  stability  analysis  has  been  performed.  Although  not  directly 
related  to  the  subject  of  this  grant,  the  first  steps  have  been  carried  out  in 
rH  ■  to  develop  a  new  decoupled  iterative  method  for  the  solution  of  the  nonlinear 
systems  arising  from  the  discretization  of  the  basic  semiconductor  equations. 

The  underlying  concept  is  to  extend  Gummel's  method  in  the  framework  of 
approximate  Newton  methods  using  singular  perturbation  ideas  to  approximate  the 
Jacobian.  In  JH  ]  numerical  experiments  have  been  performed  with  the  simplest 
version  of  this  concept.  They  show  a  significant  improvement  over  the  conver¬ 
gence  properties  of  Gummel's  method  for  unipolar  devices. 


f .  Conferences  Attended 

(1)  'International  Conference  on  Computational  Methods  in  Applied  Science', 
Paris,  France,  December  1985. 

Title  of  talk*  'Numerical  Methods  for  the  Tansient  Semiconductor  Device 
Equations',  (invited). 

(2)  'Conference  on  Semiconductor  Devices  and  Processes',  Swansea,  UK,  Duly  1986 

(3)  IEEE  conference  on  'Numerical  Analysis  of  Processes  and  Devices',  Santa 
Clara,  November  1986. 

Title  of  talk:  'An  Improved  Decoupling  Algorithm  for  Unipolar  Device 
Simulation ' . 

(4)  AMS  summer  seminar  on  ' Mathematical  Methods  in  VLSI',  Minneapolis,  April 
1987. 

Title  of  talk:  'The  Shape  of  Solutions  of  the  Semiconductor  Device 
Equations',  (invited). 

(5)  'International  Conference  on  Industrial  and  Applied  Mathematics',  Paris, 
France,  Dune  1987. 

Title  of  talk:  'Simulation  of  Tunneling  Effects  in  Semiconductors'. 


5.  Papers  Published  Under  This  Grant 


"An  Asymptotic  Analysis  of  a  Transient  P-N  Junction  Model",  SIAM  J.  App. 
Math.  Vol.  147  Nr.  3,  pp.  624-642  (1987). 

"A  Singular  Perturbation  Analysis  for  An  Asymmetric  Diode"  to  apear  in  IMA 
J.  Appl.  Math.,  jointly  with  P.  Markovich  and  C.  Schmeiser. 

"A  Singular  Perturbation  Analysis  for  the  Transient  Semiconductor  Device 
Equations  in  One  Space  Dimension",  to  appear  in  IMA  J.  Appl.  Math. 

"The  Shape  of  Solutions  to  the  3asic  Semiconductor  Device  Equations",  to 
appear  in  "Lectures  in  Applied  Mathematics." 

"Stability  for  the  Transient  Semiconductor  Device  Equations"  Zeitschrift 
Fuer  Angewandte  Mahtematik  und  Mechanik  (ZAMM)  67,  Nr.  7,  319-332,  (1987), 
jointly  with  P.  Markovich. 

"Difference  Methods  for  Transient  Device  Simulation",  in  Computational 
Methods  in  Applied  Science"  R.  Glowinski  ed.  (1986). 

"A  Second  Order  Difference  Method  for  Transient  Semiconductor  Device 
Simulation,"  to  appear  in  Trans.  Com.  Sim. 

"An  Improved  Gummel  Algorithm  for  the  Semiconductor  Device  Equations,  to 
appear  in  IEEE  Trans,  on  C.AD  ,  jointly  with  D.  Schmeiser. 
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